CPT-P036-2009 
ITFA-12-2009 



Spatial 't Hooft loop in hot SUSY theories at weak 

coupling 



Chris P. Korthals Altes 

Centre Physique Theorique au CNRS 
Case 907, Campus de Luminy 
F-13288 Marseille, France, and 

Institute Theoretical Physics 
Valkeniersstraat 65 
1018 XE, Amsterdam, the Netherlands 
altes@cpt . univ-mrs.fr 



Abstract 

The spatial 't Hooft loop measures the colour electric flux in SU(N)/Z(N) gauge theory. 
It is a closed loop of Dirac Z(N) flux and has strength k=l, 2,.., N-l. It is analyzed 
for generic k and small gauge coupling in the high temperature phase of N =1, 2 and 4 
SUSY theory up and including cubic order. In one loop order no qualitative difference 
with gluodynamics shows. However the two loop order shows a logarithmic divergence 
appear in the centre of the wall. This is because the gluinos become bosonic in the 
centre and acquire long wavelength excitations. We discuss the cure for this divergence, 
due to a gluino zero mode. The physics of Casimir scaling is explained. The cubic order 
can be obtained from gluodynamics by simple rescalings of the Debye masses and the 
Casimir parameter k/N c . 



1 Introduction 



The physics at RHIC is that of a strongly interacting plasma. In this regime QCD can 
only be studied on the lattice like the smallness of the shear viscosity over entropy ratio 
(for a review see ref. [1]). Many questions are still beyond the reach of simulations, 
like the fast rate in going to equilibrium. This is one of the reasons to contemplate 
the AdS/CFT connection and its non-equilibrium ramifications [6]. The CFT is four 
dimensional high temperature J\f = 4 SUSY, and can be studied in strong coupling by 
the calculating in the five dimensional bulk theory in weak coupling. Another reason lies 
in the Z(N) reduced effective actions [14] and their plethora of undetermined parameters 
in the the case of gluodynamics: maybe their SUSY counterparts have more to say about 
his. 

The thermodynamics of N=4 SUSY has been analysed about ten years ago [10] [11] [12] . 
These authors obtained weak coupling results for the free energy and the Debye mass, in 
order to see how they connected to strong coupling results from the AdS/CFT approach. 

Recently other observables have been studied, like the mean colour-electric flux in 
the plasma measured by the spatial 't Hooft loop [16]. 

Their results for small coupling are one-loop precision. To this order nothing qual- 
itative distinguishes the N = 4 case from gluodynamics. In this note we will compute 
the loop to the next order. This involves two-loop calculations which are remarkably 
simple in gluodynamics [8] and stay so in SUSY. However the results become in this 
order qualitatively different from gluodynamics. To two loop precision this effect is sig- 
naled by a divergence in the tension. How to cure this divergence is discussed in section 
6. This divergence is due to a gluino zero mode [15]. 

When doing the calculation it is expedient to use the multidimensional formulation 
for SUSY theories [5]. Not only N = 4, but also Af = 2, 1 cases are covered to the order 
we compute by a single parameter 6, the dimension used. 

Our results may then be contrasted with those of strong coupling in ref. [16] in 
AT = 4. 

In section 2 the problem is set up and the results for the effective potential are stated. 
In section 3 some technicalities of our perturbative method to obtain the effective action 
are given, in the next section 4 the physical reason for Casimir scaling is belaboured. 
Both of the latter sections can be skipped by the reader who is imainly interested in the 
results. The minimization of the effective potential to obtain the tension associated to 
the 't Hooft loop is done in section 5. The results are discussed in the next section and 
we wind up with conclusions. Some appendices are included to support claims in the 
main text. 

2 The set-up and the results 

Consider the plasma of a four dimensional SU(N) gauge theory with only adjoint coloured 
matter, like gluodynamics or N=l,2 or 4 SUSY. Such a theory admits SU(N) gauge 
transformations with a discontinuity in the center Z(N) of SU(N). In particular one 
can imagine a closed macroscopic spatial loop L of Dirac magnetic flux [2]. Such a 
loop is given by any static gauge transformation T4(L) that picks up a discontinuity 
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Zk = exp(ikj^), whenever we circumnavigate the loop, say clockwise. In fig.2(a) we 
show such a situation, with the loop in the (y,z) plane. 

The claim is that the quantum average of such a loop measures the mean colour 
electric flux through the loop. The quantum average is of course the Boltzmann trace, 
up to a normalizing factor: 

(V k (L)) = Tr Vfc(L) exp(-H/T). (1) 

The trace is over all the colourless, i.e. gauge invariant states in our theory. This is 
the quantity we will compute to two loop order for SUSY theories. Like its homologue 
in gluodynamics [8] it has a rapid exponential fall-off, in fact with the area A(L) of the 
loop: 

(V k (L)) = eM-QkA(L)). (2) 

The quantity Q k that controls this fall-off is called with a remarkable misnomer the 
"tension", or "dual string tension". 

Before we delve into formalism a few words on what we can expect on the basis of 
simple arguments. The spatial t Hooft loop is the SU(N) version of the Dirac string in 
quantum-electrodynamics . 

The Dirac string is engineered such that it is undetectable by quantum particles: 
the wavefunction of any particle in the theory, boson or fermion, is unaffected by the 
presence of the loop. The strength of the string is then equal to the total magnetic flux 
coming out of the monopole at its end. If it has two endpoints it describes a monopole- 
antimonopole pair. The strength of the string is the number of times the phase of the 
gauge transformation winds around the gauge group U(l). For the 't Hooft loop the 
gauge group is SU{N)/Z{N) 

The spatial 't Hooft loop V(L) is a macroscopic Dirac string L closed on itself. 
From its unobservability it would naively seem to be an academic object. However this 
needs not to be true at low temperatures, where the macroscopic quantum-mechanical 
coherence of a condensate can create Abrikosov flux tubes. Under such circumstances 1 
the spatial loop, or rather the state created by it from the ground state, will acquire an 
energy proportional to its perimeter. 

There is a less well known role played by the spatial 't Hooft loop. Mathematically, 
our closed Dirac string can rewritten through the use of Gauss' law as generator of gauge 
transformations as an electric flux loop [3] (see also next subsection) : 

V(L) = exp(i4vr / dS.TrEY k /g). (3) 

JS(L) 

Y/j is a suitably chosen colour matrix in the fundamental representation , the hypercharge 
generating the Z(N) transformation (exp(i2-7rY£;) = exp(i/c(27r/iV))l), and g is the gauge 
coupling. E is the electric field strength in the fundamental representation. Since gluons 
are in the adjoint representation of SU (N)/Z(N) their hypercharge is ±1 or 0, whatever 
the value of k. What does depend on k are the multiplicities. 

Once more one would say off-hand that no quantum particle in the theory will 
detect the area. However one might argue that at very high temperatures the thermally 

1 Not occurring in the theories we will discuss here. 
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screened particles are getting close to behaving as classical particles, escaping therefore 
the Dirac veto. This statement is correct in the limit of small coupling, where the 
classical Debye screening length Ad ~ is very large compared to the de Broglie wave 
length Xb ~ ^. 




Figure 1: Flux loop (a), Debye slab of thickness Ad with particle emitting its flux (b), 
behaviour of flux as function of distance from the loop (red) in (c) , and simplified (black) 

In that situation one can consider the slab (see fig. (1 (b)) to be a thin slab of 
width Ad, where independent particles contribute to the flux through the loop. This is 
shown in fig. (1(c).) It is then not hard to see [18] that one particle will contribute -1 
to the flux (because only half of its flux shines through the loop, giving exp(i7r)), two 
independent particles (— l) 2 , etc.. 

The Poisson distribution for a given species of screened gluons is that of Boltzmann 
particles in the fiducial volume -with on average I = tiAdA(L), n the density of the 
species-. It reads: 

P(l) = lexp(-l) (4) 
and so the flux average for this single species is: 

oo 

^P(l)(-) 1 =exp(-X D n(T)A(L)). (5) 

1=0 

Taking all species of gluons into account one has to determine how many gluons 
in the adjoint are having a non-zero charge AdYk ■ This turns out to be 2k(N — k). 
Because of the independence of the particles the exponent for one species is multiplied 
by this multiplicity: 

Qk A(L) = 2k(N - k)X D n(T)A(L) . (6) 

Qk is the quantity controlling the area-off defined in eq. (2), called the tension. It has 
nothing to do with a tension. The density n is the same everywhere inside and outside 
the minimal area of the loop. It just tells us how the fall-off of the mean flux is controlled 
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by the type of particles we have in the theory. The Casimir scaling is a signature for the 
adjoint representation. Evidently, in the confined phase we have no area law anymore 
as the colour electric flux is enclosed inside the hadrons. 

Coming back to gluodynamics, in leading order in perturbation theory this result 
(6) was qualitatively reproduced [8], with the Casimir scaling factor k(N — k) and the 
factor Xr, ~ 1/g. And indeed on the lattice [17] Casimir scaling is valid, not only at 
very high T, but surprisingly almost down to the critical temperature! Of course the 
simple coupling constant behaviour in our leading order result (6) is only valid at very 
high temperature. Lowering T will increase the coupling. For gluodynamics the two 
loop results are known [8] [18], and so are the cubic terms [19]. The cubic results are 
numerically small, but do violate Casimir scaling. Everything we know is consistent 
with the hypothesis that hard physics (on the scale T) respects Casimir scaling and that 
soft physics (on the scale gT) does not. 

For SUSY theories the gluons are accompanied by scalars and gluinos, all of them in 
the adjoint of SU(N). So we expect Casimir scaling according to the arguments above. 

In this paper we will push the perturbation series to two loop order and will argue 
that the next, cubic, order is just a suitable rescaling of results known from gluodynam- 
ics. 

In next to leading order we will find again Casimir scaling. This is not surprising 
as we argue in section 4. But to three loop order one needs a dynamical mechanism to 
explain cancellations in the planar diagram sector. 

First we want to understand how the electric flux formula (3) comes about. 

2.1 What has the spatial loop to do with mean colour electric flux? 

First let us produce an explicit version of such a singular gauge transform. 

A set of Lie-algebra matrices Yj, that generates the center through z^l = exp(i27rYfc) 
is [18]: 

Y k = 1 diag( fc, k, . . , k , k - N, k - N, . . . , k - N) . (7) 

N—k times k times 

The k are integers. 

The Gauss operator is written as: 

G a = -E\Df. (8) 

From the covariant derivative, D = d + ig[A„ we have 2 : 

G a = -E a .V + gf abc A\Et (9) 
93o = -9fabcA\Et. (10) 

The last equation defines the colour charge density gj$. 
Our singular gauge transformation is written as 

exp(iw a A a ). (11) 
2 Our notation is in terms of the SU(N) Lie algebra A = A a \ a , TrA a At = \8 a b, [A a At] = if abc \ c 
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Then we take td a \ a = a(x)Y k with a(x) equal to half the solid angle defined between 
x and L. 

The 't Hooft loop is defined as the operator version of (11): 

V k (L) = exp(i J dxG a uj a ). (12) 
We integrate inside a large sphere including the loop L the operator density V.(E a u; a ) 

( dVV.(TrEuj) = surface term on S*^ - / dS.TrEY k 2ir (13) 
J JS(L) 

The surface term is negligeable because of the fall-off of the integrand. In the second 
term the explicit form of uj is used. 
The l.h.s. is the sum of two terms: 

J dVV.(TrEto) = J dVTr(V.E)uj) + J dV.{Tr{E.V)u) 

= J dVTr(V.E - gj )uj + j dVTr(E.V + gj )uj 

We have added and subtracted the colour charge density to produce Gauss' operator 
(10). So the first term cancels on physical states. The second term equals, because of 
the normalization of the A matrices: 

J dVTr(E.V + gj p = ^J dV(E a .V + g%)u, a ). (14) 

It is the latter integral which is, up to a sign, the gauge transformation (10) and hence, 
combining eq. (13) and (14) and exponentiating we get: 

exp(— i [ dS.TrEY k ) = exp(z- / dV(E a .V + ffj >°)). (15) 
9 JS(L) g J 

in the physical Hilbert space. The r.h.s. is the gauge transformation in Hilbert space 
of exp(iuj(x)) with a discontinuity Zk in the centergroup on the surface. This was the 
definition of the 't Hooft loop. On the left hand side we see its representation as an 
electric flux operator. 

Once its role as electric colour flux operator is understood, it is intuitively clear that 
the average (1) will be an area law in a plasma [18]: 

(V k (L)) = exp(- a (T)Area(L)). (16) 

The quantity that controls the fall-off is called the electric tension, or "dual tension". 

Now we return to its quantum average (1) and discuss two other guises for this 
average. 
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2.2 Two other guises for the spatial 't Hooft loop average 

The second guise is useful for calculations on the lattice or in perturbation theory and 
is the Euclidean path integral version. The gauge transformation is replaced by the 
Polyakov loop in the Euclidean time direction r: 

,1/T 

P(A )(x)=V(exp(i drA (x,r))). (17) 
Jo 

Since it transforms under a gauge transformation as an adjoint: 

P(A$) = n(x, 0)P(A )tt\x, l/T), (18) 

the N-l phases of the diagonalized loop are invariant under periodic transforms, and 
shift under periodic modulo centergroup transformation [4]. 

So the gauge invariant content of the loop can be stored in the trace of N-l powers 
of the loop: 

P^(A ) = TrP(A ) k , k=l, ,N-1. (19) 

Imagine the Polyakov loop along a straight line in the x-direction. The straight line 
pierces the area inside the loop L. in the (y,z) plane. As long as we stay away from 
the boundary the Polyakov loop will have the same profile on any such line, and will 
only depend on z. When the Polyakov loop crosses the surface, where the loop operator 
Vfc(L) has its discontinuity, it will pick up the discontinuity because of (18). 

So what we need is the probability for a given Polyakov loop profile C(z) to appear. 
C(z) is a diagonal traceless real NxN matrix, obeying 

pW(C) = Tr exp(ikC/T). (20) 

Then we maximize this probability with the boundary conditions at infinity P(Aq)(x = 
±oo) = 1 and by the jump at the surface. 

The probability is determined by the Euclidean path integral with the constraint: 

C = U k)X 5(P( k )(A ) - Trexp(ikC/T)). (21) 

The overlined Polyakov loop is a short hand for the average over y and z coordinates. 
For this to make sense we put the system in a finite but large transverse volume V tr of 
the y and z coordinates. 

The constraint is inserted in the path integral. We then get, up to the usual partition 
function for the probability that a given profile in the x direction develops: 

J DADSDF Cexp(-S) = exp(-V tr U{C) - bulk free energy). (22) 

We integrate over gauge fields A, scalar fields S and fermionic fields F. 

U(C) is the resulting effective action. The effective action can be written in terms 
of a kinetic term K proportional to Tr(<9 2 C) 2 and an effective potential V(C): 

U(C) = J dx(^KTr(d z Cf + V(C)). (23) 
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Figure 2: Left panel: the spatial 't Hooft loop in the (y,z) plane at x=0. The arrow on 
the upper part indicates the clockwise circumnavigation defining the jump, eq. (12). On 
the right the same loop but now with the z-axis interpreted as time axis. The (x,y) plane 
contains a one-dimensional wall at x=0, and a microscopically small, periodic spatial 
r-direction along which the Polyakov loop profile P(C) = exp(i2irTa(x)Yi t ) is winding, 
with the phase a of the profile as function of x shown in panel (b). Fermionic fields are 
supposed to be anti-periodic in the r direction. 

The (y,z) plane at x=0, where the Polyakov loop jumps, is called a domain wall. This 
is a little confusing, because we think of a domain wall as a soliton, living sub specie 
aeternitatis. With only a reinterpretation of the Euclidean coordinate axes we can obtain 
such a picture, at the price of having either y or z axis as the infinite time axis. 

Hence the third guise is that of a one dimensional defect, that lives in the two 
dimensional (x,y) world with a small periodic space dimension labeled r (see fig. (2)(b)) 
. The Euclidean time is now the z-axis, so the temperature is zero. The defect is given 
by the jump the Polyakov loop looping the small periodic space dimension r feels when 
crossing the line defined by x=0 in the (x,y) plane. Thus the defect is once more a 
magnetic flux. The energy per unit length of the magnetic flux is again Qk(T). It is 
created at say z = —L z /2 and lives on in the z direction during the time L z , that we 
let go to infinity [25]. Calculationally there is no difference with the second guise, the 
Euclidean path integral (22). 

Its usefulness lies amongst more in an alternative but simpler and elegant under- 
standing of the nature of the divergencies we will meet in the kinetic part of the effective 
action. 

2.3 Results for the effective action for M = 1, 2 and 4 SUSY 

For our purpose we follow Vasquez-Mozo [13] and compute with the M = 1 SUSY action 
formulated in 5 > 4 dimensions [5]. When 5 = 10 the fermions are constrained to be 
Majorana-Weyl fermions, and we recover N=4 SUSY. When 5 = 6 the fermions are 
Weyl fermions and we recover J\f = 2 SUSY. When 8 = 4 we have by definition J\f=l 



7 



SUSY, with Majorana fermions. So generically the action we use reads: 

S = J d 4 x(^TrF* q + TrXTpDpXj p,q = 0, 1, ...,S- 1.. (24) 

The fields in our action are the bosonic ones A p , p=0,1...9, appearing in the field strength 
F pg = dpA q — d q Ap + i[A p , A q \. Both fermions and bosons are matrix- valued (in the Lie- 
algebra of SU(N), and the covariant derivative reads D p = d p + i[A p ,. 

The fermions A are spinors with 8 — 2 independent components because of the con- 
straints. Their dimension (and that of the Dirac matrices T p ,p = 0, 1, ...5 — 1) is 2 5 / 2 . 
All fields only vary with xo = r, x\ = x, X2 = y, and X3 = z, the remaining dimensions 
serve as bookkeeping device for the flavours (SU(4) for 5 = 10, and U{2) for 5 = 6). 
In bosonic computations the 5-dimensional unit matrix 1& often appears, whereas in 
their fermionic counterparts the projector matrix If = A x A (22 dimensional) with 
Trlj = 5 — 2 is present. For details see Appendic C. 

2.3.1 Some basic facts and notation 

We use throughout dimensional regularization with d = 4 — 2e, e — ► 0. 

The Feynman rules are the Euclidean finite temperature rules, with the Matsubara 
frequency Iq = 2irnT, n integer for bosons, and Iq + ttT for fermions. The constraint 
(21) will introduce a background Aq = C + gQo- Therefore the covariant derivative 
Do = do + ig[C, will remain in the in the propagators. Hence it is advantageous to go 
into the Cartan basis for all the quantum fluctuations gQ p = A p , for p > 1 and for the 
quantum fluctuations A. The Cartan basis consists of the off-diagonal matrices X lJ with 
the (ij) entry equal to l/s/2 and other entries zero. They form an eigenbasis for C: 

[C,X j ] = (C ii -C jj )W. (25) 

The potential is known, on general grounds [9], to have absolute minima at the 
locations where exp(iC/T) takes value in the centergroup. 

More specifically [18] the potential valleys of the lowest order (one-loop) potential 
V\ relating the origin to the nearest minimum with Z(N) value Zj~ are described by the 
matrices (see (7): 

C = 27rTq(x)Y k . (26) 

The scalar function q(x) is to be determined as function of x from the minimization of 
U(C). It enters the Matsubara frequencies 3 as: 

lo + 2irTq for bosons and (27) 
/ + ttT + 2-KTq for fermions. (28) 

Thus, inside the wall the role of bosons and fermions starts to get interchanged. 
In particular in the center of the wall-where q = 1/2- the fermions develop long range 
behaviour in the two dimensions left, just as the bosons do outside the wall at q=0 and 
1, but now in three dimensions. We will come back to this phenomenon in section 6. 

3 Not only in the propagators but also in the vertices! 
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The summation/integration measure is written as 



^/)" T ?/(2vr)' i -V 2£ ' ^ 

where the suffix b or f refers to whether we put the bosonic or fermionic Matsubara 
frequency. The scale \x starts to play a role in the one loop kinetic term. 

We fix the notation for the sum/integrals we meet in the one and two loop. They 
are proportional to Bernoulli polynomials and read (I 2 = (Iq + lnTq) 2 + r): 

4 log$) = B 4 (q) (30) 
(l + 2irTq) 



i 



1(b) In 



BM (31) 
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i 



4 4 = Hi) (32) 
(Z + 2vrTg) 



1(6) 



Bi(q). (33) 



The right hand sides are represented by Bernoulli polynomials B r ( defined in Appendix 
B) on the interval — 1 < q < 1. If the trace is fermionic, the argument of the B- 
polynomials is shifted by \. The integrals on the left hand side are periodic mod 1 in 
the argument q. 

Obviously one recovers by differentiation B% from etc., up to some overall factor. 
As they are periodic "polynomials" they are not analytic at q = 0, 1 

Since the breaking of supersymmetry resides in the shift from bosonic to fermionic 
traces, the non-renormalization of supersymmetric potentials makes differences like 

D r (q)=B r (q)-B r (q+^) (34) 

occur frequently. The number of bosonic (or fermionic) degrees of freedom appears 
naturally as well. We write 

A = (35) 
equal to 1, 3 and 4 for respectively M = 1, 2, and 4. 



2.3.2 Results for the effective potential 

From refs [10] [13] [12] we quote the one and two loop free energy: 

Fi = A(N 2 - 1)Z) 4 (0) (36) 
F 2 = (A~g) 2 (N 2 -l)D 2 (0) 2 , (37) 

with g 2 the 't Hooft coupling. 

The dimension of the adjoint representation appears because each member of a given 
multiplet has the same weight in the free energy. 
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(al) 



(bl) 



(cl) 




(a 2) 



(b2) 



(c2) 




(a 3) 




Figure 3: One and two loop contributions to the effective potential. Continuous lines 
are bosons, broken lines are fermions and dotted lines are ghosts. In (a) the bosonic 
contributions, in (b) the fermionic contributions and (c) the ghost contribution arc 
shown. In (d) the Polyakov loop (fat circle), renormalized by one gluon exchange, is 
inserted into the sum of the loops (al), (bl) and (c), given by the shaded blob. 



The screening masses of the gauge fields are from the same references: 

A 



m 2 D 



-g 2 T 2 . 



(38) 



The potential terms of the effective action read on the interval < q < i: 



V 
Vi 

Vo 



V!+fV 2 

2k(N - k)AD 4 (q) 



k(N - k)A 



A D 2 (q) + D 2 (0) - D 2 (0) 2 + 4B 1 (q)D 3 (q) 



(39) 
(40) 



The contributions of the individual diagrams are given in Appendix D. 

Note the factors k(N — k). They replace the factor iV 2 — 1 in the free energy and 
are referred to as "Casimir scaling". To one loop order they appear as a consequence of 
simple counting: for a given charge there are k(N-k) excitations in the adjoint with 
charge 1 and an equal number with charge -1. As the B4 (±q) take the same value we 
have Casimir scaling. 
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(a) 



(b) 



Figure 4: The contributions to the one loop effective kinetic term. The wiggly lines are 
q lines. The shaded blob in (a) is the sum of one loop boson, fermion and ghost graphs. 
In panel (b) the renormalization of the Polyakov loop is shown. 

For the two loop potential this counting is not enough as there are interactions. In 
section 4 the double line notation is used to show how Casimir scaling is nethertheless 
realized. 

The factors A = count the number of 5 dimensional loops. The last term 
of V2 is the exception, for a simple reason. This latter term is due to the one loop 
renormalization of the Polyakov loop (see fig. (3 (d)), which only involves the gauge 
sector and is proportional to By. 



5q = - 



~9 2 



(4vr) : 



-(3-0£i(?) 



(41) 



Inserted into the result for Vi(q), V±(q + Sq), it gives a result linear in and B3. We 
show in more detail in section 3 how these insertions come about. 

To obtain the old result for gluodynamics, put 5 = 4 and drop the fermionic contri- 
butions to find [8]: 



Vy + g 2 V 2 = k(N - fc)^TY(l - q) 



4vr 2 



-"4 2/ 



(42) 



2.3.3 One loop kinetic term 

We now turn to graph (a) in fig. 4. It is obtained by computing the determinant of 
boson ghost fields in (69) and expanding the determinant in terms of derivatives of the 
background. We suppose the colour degrees of freedom have been diagonalized in the 
Cartan basis Q 13 , so that the background B below stands for B %3 = Ba — Bjj. 

So we have to take into account that the x-dependence in the background does not 
commute with the momentum operator P x = jd x in the d'Alembertian. The B field 
will depend on the operator X, B_ = B(X), that obeys: 

= — i and 

X\x) = x\x). (43) 
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We do this in Feynman gauge £ = 1. This gives for the bosonic part of the effective 
kinetic term ( for the fluctuation matrix see (69): 

K b = Tr fe ^ log ( - D 2 (Bl h - 2F 0a; lox) - Tr fe log ( - D 2 (B_)j . (44) 

The first term is due to vector bosons and scalars, the second logarithm due to the 
ghost. The magnetic moment term of the gluonic component is proportional to Fq x = B' 
and hence O(g) and expanded out. The scalar component has no such term, as reflected 
in the matrix: 

l Ox = Diag(l,-l,0,...,0). (45) 

So Kb becomes: 

Kb = S -^Tr b \og(-D\B)) - ±F 2 x 4Tr b { _ D l {B)f ) • (46) 

We do the same for the fermionic determinant. The result there is different in that 
the spin terms are present for all fermionic components 4 . 
One gets accordingly: 

K f = -Tr f ^log(-D 2 (B)1 { -F 0x r 0x lf) 



Tr f \og{-D{BY - y KTr f { _ D{m ) . (47) 

The spin matrix T 0x = ^[T , T x ]. 

The trace is mixed in x and conjugate momentum, instead of the pure momentum 
trace J^in section 2: 

r d d ~ 2 l r 

TrO / J^Y^Te j dx(x\0\x). (48) 

The suffixes b or f refer to treating the Matsubara frequency lo = 2irnT as bosonic or 
fermionic, by adding in the latter case ttT. 
The total effective kinetic term is now: 

K = K b + K f (49) 

and from (46) and (47) follows: 



K = ^(Tr s log(-D 2 (B))-T r/ ;o 9 (-D 2 (B))) 

+ kM- iTn rikw + ^ Tr 'rikw) m 

The short distance behaviour of the total effective kinetic term is independent of 
the background and whether the trace is bosonic or fermionic . The short distance 

4 Easily checked by going into an explicit representation for the Go and G x matrices. See also 
Appendix C. 
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contributions of the first two terms in Kb and Kf 5 do therefore cancel in the effective 
kinetic term, and we conclude the short distance behaviour of the kinetic part of the 
SUSY action is fixed entirely by the difference of the magnetic moment terms of bosons 
(—4) and fermions (^f^): 

1 (5 2 1 

Kshartdistance ~ ^ x ((-4 + ~^-)(- + 0{e)). (51) 

This difference vanishes for 5 = 10, as it should. If S is smaller than 10 we retrieve 
the asymptotically free SUSY theories. For S larger than 10 asymptotic freedom is lost. 

Finally, from (118) in Appendix A and (50) the kinetic term equals in d —> 4 and 5 
dimensions: 

K = 

9 Z 

g 2 r(^) fT\- 2e (f^ ... 2e5-2\^ 1 



2e<S-2\ v ^ 1 1 g 



x2 



(52) 



|n + 1/2 + g| 5 - rf J (4vr) 

The very last 0(g 2 ) term comes from diagram (b) in fig. (4) and results from substituting 
the renormalization 5g = —j£^yz(3 — £)Bi(q) ( see eq. (72)) into the kinetic term q' 2 in 

Feynman gauge £ = 1 6 . Details can be found in section 3. The one loop /3-function 
coefficients bb respectively bf are the contributions of the bosons respectively fermions: 

b b (S) = + ^ (53) 

bf(S) = (54) 

The total /3-function is then bo(6) = bb{5) + bf(5) = | — 5, asymptotically free for 
5 < 10. 

Substitute [24]: 

in (52). 

In terms of the renormalized coupling g(T) 

(56) 



2 2 (T) 5 : 
and (5, we obtain for the kinetic term: 

(2^' Or)) 2 / /^(T)^ 2 



- — log 7T — 2 log ( — 



+ + 5)) -WW(l/2) -■»})• (57) 



5 Working out their contributions proportional to Fq x = B' 2 is straightforward and is done in Ap- 
pendix A 

6 Including this renormalization the £ dependence drops out from the kinetic term as discussed in 
ref. [8] and [20] 
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The representation in terms of the logarithmic derivative of the T function ip (see 
[24]) is valid for < q < \. The reader will find the details in section 4 and Appendix 
A. 

3 Steepest descent of the constrained path integral 

This section will deal with the constraint (21) in the Euclidean path integral. The main 
point is the systematic approach, which is of import for higher loop calculations. For 
readers who are only interested in our two loop results this section can be skipped. 
We reproduce for discussion's sake (22): 

J DADSDF Cexp(-S) = exp(-V tr U(C) - bulk free energy) (58) 

and the constraint C: 

C = Jl ltX 6(P(i)(Ao) - Trexp(iZC/T)). (59) 

The steepest descent expansion of the Euclidean path integral is straightforward, 
and in principle there is nothing new except the handling of the constraint C, in (59). 
This new feature is discussed in this section, in particular the diagrams in figs. (3 (d)) 
and (4 (b)) in the second subsection. 

3.1 Generalities on expanding the constraint 

The most obvious way of dealing with the constraints ( which are only depending on 
the space variable x) is to Fourier analyze them with one-dimensional x dependent fields 
7®(x) where 1 runs through the integers from 1 to N-l. Since the constraints are gauge 
invariant, so are the constraint field 7. We then get an action Scan, with these constraint 
fields included, 

S C on = i E / <W° (a) (pW(A q )(x) - exp(ikC(x)) + S. (60) 

The constrained action appears in units of the squared coupling g 2 in the path integral: 

/ DADSDFDjexp(-^). (61) 
J g z 

We start the steepest descent around the field C(x) by introducing background fields 
B(x) and jc\x) in the potential and the constraint fields: 

A = B{x)+gQ , (62) 
7 (/) (x) = 7 P(z)+37? ) (*). (63) 

B(x) is diagonal, traceless of size N x N. The subscript q refers to the quantum 
fields that fluctuate around the gauge invariant sources 7P (x) and Qo fluctuates in four 
dimensions, around the background field B. The lattter is not gauge invariant. 
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The expansion in the quantum fields in S CO n yields equations of motion that fix 7P 
and B in terms of C, and define a minimum of S con . They read respectively: 

(p®(B)-exp(ilC/T)J = 

i^lP^(B)+d x 2 B d = 0. (64) 

Derivatives with Q d {x) are indicated with subscripts , d. Only diagonal and over lined 
Q's are relevant, because B is diagonal and both B and 7W do only depend on x. The 
matrix (C) = IP® (B) brings us from the d-basis, numbering the diagonal elements 
of the Lie-algebra, to the 1-basis, numbering the powers of the unitary Polyakov loop. 
It reflects the non-linearity of the constraint. Surprisingly, in amplitudes it will cancel, 
as we will see in the next subsection. 

So B=C from the first equation, and the second equation says the d component of 
i^jp equals the d-component of the double gradient of B (or C). From the first linear 
equation follows that only 

^Tr^C) 2 (65) 

remains as classical term. So there is no classical potential term. This has an important 
consequence: the potential comes entirely from quantum corrections, which are of O(g ). 
As a consequence the gradient of the C field must be soft: 

d x C(x) ~ g. (66) 

So the second linear equation in (64) tells us that the gauge-invariant source is soft 
as well; the 7P are 0(g 2 ). 

Turning to the quadratic terms in S con we have in linear £ gauge the following terms: 

Q M ( - D 2 (C)S^ + (1 - ^)D^D V - 2i[F^Q v + ghost terms (67) 

+ ij l q k) kP^(C)Q d (68) 
+ i 7 VpV\C) >ab Q a Q b (69) 

We write in the sequel QW = kP®(C)Q° d . 

We have kept in these terms all C dependent terms, although their derivatives are 
soft, like the last term with the gauge invariant source. 

The first term is the familiar inverse background propagator, the second one, after 
integration over the constraint fields, gives us delta function constraints on the Q®. 

(k) 

The third one, because of the softness of the source 7c , will give a one loop correc- 
tion to the classical term (d z C) 2 as we will see in the next subsection and is entirely due 
to the constraint. It introduces a renormalization of the background field in the kinetic 
term. 

The new feature entering in cubic and higher order terms is the appearance of terms 
proportional to ijg \ e.g. ij^ P^QaQb- 
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If we want to integrate over the fluctuations 7^ to reinstaure the constraint on the 
quantum fields Qq we can replace the factors; 

i 7 W exp(i 7 (')Q(0) ^ _JL exp(n«Q(0) (70 ) 

Integrating now over the j q we get the constraints on the Qo- The partial derivative 
— ^Tjy acts through partial integration on the vertices of the usual Yang-Mills action with 

gauge fixing and ghosts and creates new, "inserted" vertices. Such vertices are shown 
as fat dots in fig. (3 (d)). 

3.2 Insertions to two loop order 

Let us now look at what these insertions of the renormalization of the Polyakov loop 
imply for the two loop calculation. 

We start this with a useful identity [20] for this renormalization, and the two cor- 
rections to order g 2 that follow immediately from this identity. 

The Polyakov loop is shown in fig. (5) with one gluon going across. This is the lowest 
order g 2 . On the righthand side are shown the ways the gluon can go across. Ignore the 
partial derivative with Qq^ for the moment. If one adds all the possible exchanges the 
result is : 

(P%(C)QlQl) = t«\C) 4 5 2 C d , with: (71) 
faCu - 5 2 Cjj =- -^B^CiW-Z). (72) 

The first identity causes the simplicity of the one-loop insertions, because all the 
non-linearities of the constraint are sitting in the transition matrix t d l \c). In fact, as 
promised, it will disappear from the insertion graphs in fig. (3) and (4), once we look 
at the graph whole. 

For the insertion graph in the potential term the partial derivative acts on one of 
the gauge field vertices. But since we sum over 1 the transition matrix in Q$ cancels 
with the one in (72) and we are left with 5 2 C d -S^ for the insertion graph. 

To find the insertion graph for the kinetic term apply our identity to the third term 
in (69) and combine the result with the second equation in (64): 

i4 l) (P {l) (C), ab Q a Q° b ) = -(3 - Od 2 x C d .5 2 C d = -(3 - ^CM^B^dffZ) 



= $-Z)T±[d x B,[\j\ji]]d x B x (Cij) (74) 

= (3-0^QA£i(Cy) (75) 
Now use that B\ ( given in appendix B) is linear in its argument : 

Bi(B ij ) = ~(B ij -^e{B ij ). (76) 



The derivative of the first term is what we retain in the graph. The derivative of 
the second term equals 5(Cij)d x Bij. It can be ignored in the tension because of the 
equations of motion for B (see section 5) . 
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Figure 5: One gluon exchange in the 1-th power of the Wilson line. Adding up all the 
terms inside the square bracket gives eq. (72) in the main text. 

In summary: the insertion amplitudes are free of the non-linearity. They amount to 
a linear renormalization of the C matrix, as in the second equation of (72). This means 
that they drop out when determining the value of the minimum of the effective action 7 

But they are important for the gauge independence of the effective action, and hence 
for the equations of motion that determine the profile, i.e. the location of the minimum. 

4 Effective potential to two-loop order and Casimir scaling 

Casimir scaling has been proven at the two loop level in a brute force way [19] [16]. 
Although not of direct relevance to our results we think it may be useful to remedy the 
situation. 

To understand the Casimir scaling of the potential in a transparent way at the two 
loop level, we introduce the double line representation by 't Hooft, valid for finite N, 
and with a simplification for the class of models at hand. This is important because 
we will thus see that Casimir scaling is not only a consequence of counting of degrees 
of freedom, but maybe a genuine property of the hard modes in an interacting plasma. 
From the method below it will be clear that for three loops Casimir scaling necessitates 
constraints on planar graphs. Non-planar three- loop graphs are Casimir scaling. 

The double line representation is based on the fact that the adjoint representation 
can be generated by the action on an N x N Hermitean traceless matric Q of the 
fundamental representation on the right and from its contragredient partner on 
the left: 

AdQ(Q) = QQQ- 1 . (77) 

Remember from previous sections the Cartan basis in group space. This basis was 
instrumental in diagonalizing the fluctuation determinant in the background of the Wil- 
son loop. We recall its general form 8 . It is a set of N x N matrices X 1 ^ , where the labels 
i and j run both independently from 1 to N. They are given by their matrix elements: 

tftl = (78) 

7 The value of the minimum is computed as an integral over C 

8 Some of the material below was developed in conversations with Rob Pisarski 
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They form an orthogonal set with norm 1/2: 

Tr\ ij X kl = h il 5 jk . (79) 

These elements form a basis for the U(N) algebra. This is obvious from the fact that 
i and j run independently through 1 to N. The N diagonal X n are an alternative to the 
traditional diagonal basis: 

\ d = * diag(l,l, ,l,l-d,0,...,0),d = l,...,JV-l (80) 

A = -4=1- (81) 

To avoid the U(l) component, one has to change the diagonal members of the Cartan 
basis by: 

Am = A" - -^-1 (82) 

so that they become traceless. 
In terms of matrix elements: 

4 = - ^6 kl ). (83) 

Now the fluctuation fields Q carry the indices ij when we decompose the NxN matrix 
Q = Q l ^\ l K And therefore the propagator {Q l ^Q nm ) carries four indices and becomes: 

(Q ij Q nm ) = (S im 5 jn - -^5 ij 5 nm )/{p ij ) 2 . (84) 

We did not put in the Lorentz indices for notational convenience. 

In fig. (6) we put the Kronecker delta's in evidence by the lines connecting the 
indices. 

The vertices in the Cartan basis are also products of Kronecker delta's. The reason 
is that the structure constants are formed by commutators, and products read: 

X ij x im = J_ 5 ji X im . (85) 
v2 

Together with the orthogonality relation (79) this gives the result: 

2y/2f i i' lm ' nr = 5 ri 5 mn Si l — p n S rl S mi . (86) 

In the figure we put only one of the two terms. The relative minus sign is important 
when comparing planar to non-planar diagrams. 

The second term in the propagator is the "U(l) gluon" , proportional in this notation 
to the unit matrix 1. This gluon does not couple into any of the vertices, since the 
structure constants are coming from commutators. For the propagator of the gluino the 
same is valid. 

Therefore in our problem we can just as well drop the second term in the propagator. 
This gives rise to a big simplification. 
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Figure 6: Double line notation for propagators and vertices of bosons, based on (84) 
and (86) in the text. For fermions and boson-fermion vertices analogous rules are valid. 

From the double line representation for the propagators and vertices it follows imme- 
diately that the colour shift in the Matsubara frequencies is conserved. Take for example 
the three gluon vertex, and draw a circle with midpoint in the vertex. Then the shifts 
in the three outgoing legs add up to: 

C ij + C ji + c u = ( 87 ) 

Consider any one loop diagram for the free energy. We suppose it is evaluated in 
the background C = 2-rrTqYk and we recall from section (2), eq.(7), the form of 

Y k = — diag(fc, k, . . . ,k,k — N,k — N, . . . ,k — N). (88) 
N v ' v ' 

N—k times k times 

The k are integers. The rows and columns of are indexed by i,j,...If i,j,.. are giving 
one of the first N-k elements such and index is said to fall in the "k-sector" [k]. If the 
index gives one of the last elements it is said to fall into the sector [N-k]. 

Clearly, if the indices (ij) of a given double line are in the same sector the corre- 
sponding shift in p l Q = po + q((¥k)u ~ (Xk)jj) 1S zero - If they are in different sectors the 
shift is ±q. So for a one loop diagram there are 2k(N-k) possibilities of having q flowing 
in one or the opposite direction. This because there are k ways we can chose hte index 
i from the [k] sector and independently N-k ways to chose j from the [N-k] sector. The 
weight of any of those configurations is the same, ~ B^^q), and so we get that sum is 
proportional to k(N-k) which is called Casimir scaling. 

In the rest of this section we will show that any two loop diagram will display Casimir 
scaling. 

For this we recall that the well-known convenience of the double line notation comes 
from the case where q=0 (no background). The diagrams form index cycles, and every 
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cycle is counting for a factor N. This leads then to the topological classification of double 
line diagrams into planar, and non-planar diagrams. An example is shown in fig. (7). 
Clearly the right most diagram has two index cycles less than the planar ones on the 
left. 

However if q is non-zero the counting becomes a little more complicated, and the 
rest of this section is handling this problem. 




no q-line one q-circuit non-planar diagram 



Figure 7: A two loop diagram in double line notation. 

Let us turn these simple observations into a proof that any two- loop diagram with 
non zero q carries a factor k(N-k). 

Look at fig. (7). Any line (propagator) that carries a Matsubara frequency (no + q) 
is called a q-line. So the diagram falls into two classes: one where there is no q-line, or 
one where there is at least one q-line. The non-planar 2 loop diagrams have only one 
index loop (they are 0(1/ N 2 ) smaller than the planar loops with three index cycles). 
So they have no q-lines. 

The planar 2 loop diagrams can have q-lines as shown in fig. (7). 

Remember in a vertex there is conservation of colour and hence of q. So in the two 
three-vertices there is an ingoing q-line and an outgoing q-line. Hence a q-line follows 
a closed loop, in analogy with the chemical potential for fermion lines. An example of 
such a closed loop is given in fig. (7). It is indicated by the broken line in the second 
diagram on the right. 

We are now in a position to tell how many of those q-line diagrams there are. To 
have a q-line the two indices constituting the line must be in different sectors of the 
matrix. In the diagrams in fig. (7) we must have m and i in different sectors. And j 
must be in the same sector as m, and in a different sector from i. 

So we have the following situation in table (1): 

The first two lines in the table correspond to the planar diagram with the q-circuit 
as shown in figure (7). The next 4 lines to the two planar diagrams with the q-circuit 
drawn differently. The last two lines correspond to no q-lines. 

The first line has (N — k) 2 possibilities for the indices m and j, and k possibilities 
for the index i, so k(N — k) 2 in total. The second line has k 2 (N — k) possibilities and 
describes the diagram with the same q-cycle as in the first line. So the total is: 

k{N -k) 2 + k 2 (N -k) = k(N-k)N. (89) 

The factor N is absorbed by the coupling g 2 , and so Casimir scaling is verified! 

For q=0 we should find the familiar factor N(N 2 — 1) for the free energy. The 
terms just discussed contribute 3k(N — k)N. The last two lines in the table contribute 
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Table 1: Counting of the multiplicity of the two-loop planar diagram, [k] is the Yk sector 
with N-k entries k, and [N-k] is the sector with k entries N-k. i, j and m are the indices 
of the index cycles in fig. (7) 



[*] 


[N-k] 


mj 


i 


i 


mj 


3 


mi 


mi 


3 


m 


ij 


ij 


m 


ijm 


ijm 



A; 3 + (N — k) 3 . The contribution of the non-planar graph comes with a relative minus 
sign 9 and equals — (k + (N — k)), because there is only one q circuit. The sum is indeed 
the familiar free energy factor. 

Prom this analysis one can draw an immediate conclusion for the 3 loop potential: 
the non-planar contributions with only two (or less) q cycles will show Casimir scaling. 
But the planar diagrams necessitate relations between weights . 



5 Minimizing the effective action 

From the explicit expressions for kinetic and potential terms in the effective potential 
developed in the previous sections we now find the mimimum of the effective potential 
by varying the profile q: 

U = J dx{K + V). (90) 

The kinetic term is proportional to q' 2 , K = Kq' 2 . Then the minimum of U is found 
by quadrature: 

Umin = j dx(y%J - W^j +2^ dq\[kV. (91) 

The first term is non-negative definite and produces the equations of motion for the 
profile q by setting it zero. The second term is the value of the minimum. 

There are two poles in the kinetic term. One is at q=0 and is canceling with the 
corresponding zero of y/V ~ q(l + regular in q). But the pole at q = 1/2 in ip(l/2 — q) 
causes a logarithmic divergence. This pole is due to a single Matsubara frequency in the 
sum for the gluinos. In fact the denominator |n + \ + q{x)\ in eq. (52) becomes bosonic 
at the center of the wall, where q(0) = 1/2. So the gluino becomes at the very center of 
the wall a boson. 

9 Vertices in non-planar graphs get their contributions from terms in the commutators which have 
opposite sign 
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And its Matsubara frequency n = — 1 causes the divergence. We will just by fiat 
subtract it so replace Y>(l/2 — q) by 

+M-*>-Tjk? (92) 

The divergence has an interesting physical explanation, and we come back to it in 
section 6. 

To wet the appetite of the reader: the remedy is not simply to reinstall the thermal 
mass of the gluino in the subtracted pole in eq. (92). It would help to cure the divergence 
because 

,1/2 1 

d lT7^ 7 T7W\ ~ !og(m(? = 1/2))- (93) 



/ 



l/2-q + m x (q= 1/2) 

However, the thermal gluino mass turns out to be zero for q = 1/2! 

To get the tension we plug in the results for K and V from eq. (57) and (40), and 
using (92) one finds: 



g k (T) _ 2V2tt 2 ( 9-2V3) A ^r 



1_ jty{ {( ~ 2 - 92683 - ] xbb{6) 



k(N - k)T 2 15 g 

+ (3.27471...) x b f (5) + (1.96351..) x ^} + x 5 J 

For gluo dynamics we have: 



(94) 



4tt 2 T 2 

g k = k(N-k) — =x 
3V35 2 

~~2 



1 " -J^{{(-2.98505... + 1.96351..) x (-H) _ I _ 2 } + I x sj 



(95) 



The terms within small curly brackets stem from the kinetic tern folded with the 
square root of the lowest order potential. 

Numerical values for ^(1/2) = —1.96351.. [24] and for the folding of the psi-functions 
with the square root of the lowest order potential are used. 



6 Discussion 

First we justify dropping the gluino mode that gets bosonic in the center of the wall. 

To this end we recall the third guise for g k as the energy of a one dimensional soliton 
as explained at the end of section 2. The soliton is stretched along the y-direction, at 
x=0. It separates two regions in the x-y plane with the Polyakov loop having the value 
1, respectively exp(ifc 2 ^). 

The reason we can drop the mode that diverges at x=0 is that it has zero energy 
in the field of the soliton. Remember that the time is the z-direction, so zero-energy 
implies for the gluino spinor: 

r z d z x = o. (96) 

To find the zero-energy solution, write the Dirac equation for the gluino, with the 
momentum in the r direction irT — 2nTq{x):. this is the mode causing the divergence 
in i/>(l/2 — q) for q = 1/2 in the effective kinetic term, eq. (57). 
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We have only x dependence left for the zero-energy mode: 

(r x d x + T (itT - 2irTq(x))^j X(x) = (97) 
So the solution is [15], u being a constant spinor: 

Xq(x) = exp (t T x J dx'(-nT - 2irTq{x')) S ju. (98) 

To have a localized state Ao we need YqY x u = u. 

We have found that the diverging mode is a zero mode. It binds to the soliton and 
forms a state degenerate with it. Hence it is legal to drop it. 

There are several interesting aspects of this state that will be discussed in a sequel 
paper. 
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Figure 8: Behaviour of the half-integer gluino energies due the presence of a magnetic 
flux. Vertical axis is energy in units of 2nT, horizontal axis is the x-axis, orthogonal to 
the magnetic flux in the y-direction at x=0 (see fig. (2 .b). Gluinos are anti-periodic 
in the short spatial direction r, and have zero momentum in x and y direction, as in 
eq.(97). The continuity of the derivatives at the jump at x=0 is manifest. 

We mention as a curiosity the contribution from the potential to the tension: for 
gluodynamics it is from eq. (95) 5/2, and for N = 1, 2, 4 it is given by ^5. In the case 
of gluodynamics this is already true on the level of the effective potential, eq. (42), for 
any value of q as in eq. (42). In the case of SUSY it is only true for the minimum Qk , 
i.e. only after integration over q. 

The kinetic part is subject to ambiguity due to the choice of the running coupling. 
With our choice the 0(g 2 ) correction is negative and qualitatively like in gluodynamics. 
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6.1 The profile to lowest order and the fluctuation operator 

Up till now we did not look at the other piece of information the effective potential holds 
in store: the profile corresponding to the tension. 

The profile follows from the first term in eq. (91), substituting (40), (34) and (125) 
and (129). The equation of motion becomes to lowest order in terms of q and £ = mr>x, 
mo the Debye mass from eq. (38): 

q' = ^1-y (99) 

Note that the profile does not depend on k, the strength of the loop, and that 
this equation is valid for < q < 1/2. Continuity of the kinetic term and symmetry 
arguments extend it to the upper half of the interval. Clearly the profile is for small q 
exponentially decreasing as exp(— tod|x|). 

The solution for the profile is: 

q m = ^(l-tanh 2 (-3m D x/8 + x )) (100) 

with xq such that q(0) = 1/2. 

The 0(g 2 ) correction to the profile is straightforward as long as we stay away from 
its wings (small q). The reason is the pole at q=0 in the kinetic term. This pole signals 
the well-known presence of corrections to the Debye mass which are not 0(g 2 ), but 
larger, by a logarithmic amount. In ref. [19] this problem is discussed in more detail. 

A last remark concerns the 0(g 3 ) corrections. They can be obtained from their 
gluodynamic analogue, as the following argument shows. 

The second derivative of the lowest order potential V\. It reads: 

V{'(q)=m 2 D (l-4q). (101) 

This turns out to be the one loop self-mass of the time component of a stationary, 
zero spatial momentum gluon, in the background q. The gluon itself is supposed to 
be q-neutral, so is not screened by the background. Such gluons are the cause of the 
familiar, electric infra-red divergencies present in the plasma without background. 

The q-independent part of the self-mass gives upon resummation the Debye mass 
corrected gluon propagator, and hence the 00(g 3 ) correction to the free energy: 

Free energy due to soft modes ~ 1/2 log det(— V 2 + m 2 D ) ~ m|> ~ g 3 . (102) 

Taking into account the q dependence of the selfmass one computes the same cor- 
rection but now to the tension through the determinant. Since we now compute soft, 
O(gT), fluctuations, which are of the same order as the length scale in the profile the 
use of gradient expansions is not allowed in this determinant. So one has to evaluate 
the eigenvalue spectrum of the Schroedinger equation: 

tension due to soft modes ~ ^ logdet ^ — V 2 + m 2 D (l — 3/ cosh 2 (— 3mDx/<$ + xo))^J • 

(103) 

However this calculation has already been done in gluodynamics[19], where the fluctu- 
ation operator V{'(q m ) is identical, up to some trivial rescalings. 
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7 Conclusions 



In this paper the corrections to the tension of the 't Hooft loop for hot SUSY theories 
have been computed to two loop order. The cubic order is shown to be a fairly trivial 
rescaling of the same problem in gluodynamics. The calculation can be done analytically 
because the one-dimensional Schroedinger problem, eq. (103), is separable. 

What can the understanding of the spatial 't Hooft loop bring us in terms of a 
better understanding of the hot phase? After all, the most we learn from it seems to be 
that it scales according to the Casimir scaling law, k(N — k) and that this scaling law is 
characteristic for gluons, the screened quasi-particles, being in the adjoint representation. 
This comes hardly as a surprise. 

However, in J\f = 4 theory the self-duality of the theory at zero temperature may help 
to shed light on the behaviour of the spatial Wilson loop. The latter measures the mean 
flux of magnetic quasi-particles . It is well-known from lattice calculations [26] that 
also the Wilson loops show Casimir scaling. It is tempting to associate to this scaling 
magnetic quasi-particles in the adjoint representation in some magnetic group [18] [22] [21] 
and Af = 4 theory seems a setting where one may have a more quantitative grip on these 
ideas. 
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Appendix A: Determinant in mixed x and p representation 

In this appendix we elaborate on how to compute the determinant of the d'Alembertian 
with a background B, where the variation with x does not commute with the corre- 
sponding derivative. We suppose the colour degrees of freedom have been diagonalized 
in the Cartan basis Q lJ , so that the background B below stands for B tJ = Ba — Bjj. 
We start with eq. (44) in the main text and write 



for the difference between the non-commuting (underlined background B) and the com- 
muting case ( written as B). 



Trlog(- D 2 (B)) - Trlog(- D 2 (B)) 



(104) 



We write: 



B = B(X). 



(105) 



25 



The operator X and the momentum operator in the x-direction have the commutation 
relation as in eq. (43) 

[X,P x ] = -i (106) 

and they act in a Hilberspace spanned by the normalized kets \x), with X\x) = x\x). As 
in the main text we introduce the eigentime t and write the difference of d'Alembertians 
as: 



-Tr 



exp(-t(-D(B 2 )) - exp(-t(-D(Bf 



dt 
T 

The trace becomes in this language (Iq = 2nTn, n integer) 



(107) 



Tr exp(-t((l + B) 2 + P t l + P^ 

= T^J^^^Jdxixlexp^-mo + BiX^ + ll + P^lx). (108) 

The background field B(X) = B(x + (X — x)) is expanded around B(x), the field 
that commutes with P x : 

B(X) = B{x) + B'{x){X -x) + ^B"(x)(X - x) 2 + ... (109) 



We know already- from the absence of a classical potential term- that every derivative 
B'(x) = dB J^ counts as a factor g, and since we are interested in order g 2 we neglect 
the dots. 

The matrix element in (108) becomes on substitution of the expansion for B(X) ( 
we use the notation 2A = ((/ + B{x))) 2 )' and = A': 

(x\ exp(-i(/ + B(X)) 2 + l 2 r + P X 2 ))| X ) 
= <x| exp(-t(/ + B{x)f + n(X-x-^) 2 -^ + l 2 tr + P 2 ))\x) 

= (^\ex P (-t(l + B(x)) 2 + n(X) 2 -^+l 2 r + P 2 ))\^). (110) 

The dependence on x in the last expression has been translated away. It stays only 
implicitely in B(x). 

In what follows we assume it makes sense to expand the integral (108) in terms of 
A and Q. 

We note that in case ft = the matrix element becomes proportional to: 

( x \exp(-tP 2 )\x) = -^- T . (Ill) 

(4vrf)2 

This is precisely the result one would get from the l x integration, pretending the back- 
ground is constant. 

In our case however we have Q, = 0{g 2 ) and borrow a formula from [23] to take its 
presence into account: 

<4ll«xp<-.(tur» + f2))|j*) = ^«p(-i£ ( -i + raM ) (112) 
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with 

a = (113) 
sinn a; 

uj = 2Dt. (114) 

In the second term expand a and coshw for small u and substitute the result in 
(110). Terms with Q, can be transformed 

J dxexp(-t(l + B(x)) 2 )Q = J exp(-i(Z + B{x)) 2 )2A 2 t (115) 

through partial integration in x and the definition of A and $7 above eq. (110). The 
boundary terms at x = ±oo cancel with each other for B(±oo) = 0, or 2nT after 
summing over all n in Iq = 2irnT. The summation over Matsubara modes introduces 
the periodicity in B mod 2irT, hence the cancellation. 

The result of this expansion becomes for the matrix element in (114) simply: 

(A, exp( _ ((si ^ + Pm A.) = -Lyd _ l A V + ..). die) 

Looking back at eq. (108) we see that the integrations over l tr produce a power 

2-d 

of eigentime T 2 . Together with the powers of t coming from the expansion in the 

3 , 5-d 

oscillator frequency f2 above , i.e. , one arrives at a power t 2 . 

This power ensures that the term with the factor A 2 = (lo + B{x)) 2 B' 2 becomes, 
after substitution into (107), through partial integration in the eigentime t: 

A 2 = (117) 

hence our final result for (104) becomes, using (107), (116) and (117) : 
Trlog(- D 2 (B)) - Trlog(— D 2 (B)) 

We wrote q = B(x)/2irT, and e = 2 — d/2, and to complete the result we take from 
[24] for d -> 4: 

E^ = ^ (?) ^ (1 " 9) ' (119) 

The i/j function is the logarithmic derivative of the Euler T function. It has a simple 
pole at q = 0. 

In the main text ( eq. (44) and (47) ) expressions of the form: 

Tvlog(-D 2 (B) ± B') - Trlog(-D(B)) (120) 

appear. 

Expanding the B' gives us for half the sum of the terms with plus and minus sign: 
/ dX ^2~ Tl {-D{B)Y + T ^g(-D(B) 2 ) - TTlog(-D(B) 2 (121) 
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By the same technique the first term leads to eq. (118), with the factor \^- 2 ^ L replaced 

by-i 

The difference between boson respectively fermion trace is the periodicity respec- 
tively the anti-periodicity. Hence it comes into play in the argument of the ip functions. 
For bosons the argument is q, for fermions q + ^. 

The former produce a pole at q = 0, the latter a pole at q = 1/2, as discussed in the 
main text. Both these poles signal the presence of infra-red effects, but are remedied in 
different ways. 

Appendix B: Bernoulli polynomials 

They are related to the polynomials defined in (33): 

B d=4 (x) = ^ir 2 T 4 B 4 (x) (122) 

B 3 (x) = It:T 3 B 3 (x) (123) 



B 2 (x) = \t 2 B 2 {x) (124) 
M*) = ~TZ B M ( 125 ) 



and finally (e(x) is the sign function): 

B A {x) = x 2 {l-\x\) 2 (126) 
B 3 (x) = x 3 - -e(x)x 2 + -x (127) 

B 2 (x) = x 2 -e(x)x + ^ (128) 

B 1 (x) = x-^e(x) (129) 

Appendix C: Spinor projectors in 5 dimensions 

We establish the form of projectors on Major ana- Weyl spinor s in 10 dimensions, on 
Weyl spinors in 6 dimensions and Majorana spinors in 4 dimensions. This is important 
for establishing the dimensional factors in the fermion amplitudes. 
•6 = 10. The projector If = XX is defined by A = A C and TuA = A. 

• 5 = 6. The projector is defined by ^A = A 

• 5 = 4. The projector is defined by A = X T C 
In all these cases one has Trlo^l/ = 5 — 2. 

Appendix D: Individual diagrams contributing to the po- 
tential 



We compute the potential term V = V\ + g 2 V 2 in the effective action: 



U = dx 



K(q)q' 2 + V(q) 



(130) 
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We refer to fig. (3) and compute in Feynman gauge. All two-loop results below are 
given in terms og the 't Hooft coupling g 2 = g 2 N . The hatted Bernoulli polynomials 
are given in Appendix B. 

We get for the bosonic graphs (a): 

ai = 5 -2k{N - k)B A {q) 

a 2 = -^(5-l)k(N-k)^B 2 (q) 2 + 2B 2 (q)B 2 (0^j 

a 3 = - A 5{5 -l)k(N -k)(^B 2 (qf + 2B 2 (q)B 2 (^j. (131) 

For the ghost graphs (c): 

ci = -2k(N - k)B 4 (q) 

c 2 = ^k(N-k)^B 2 (q) 2 + 2B 2 (q)B 2 (0)y (132) 
For the one loop fermion and mixed boson-fermion graphs (qf = q + 1/2): 
bi = - 5 -^2k{N -k)B 4 {q f ) 

h = (^) 2 HN~k)\ l B 2 (q f ) 2 + 2B 2 (l/2)B 2 (q f ) 

- 2 (B 2 (q f )B 2 (q)+B 2 (1/2)B 2 (q) + B 2 (0)B 2 (q f )^ (133) 

Compare to eq. (40) and (34) in the main text. Adding a\,bi and c\ gives V\ . V 2 
is obtained by adding a 2 ,a^, b 2 and c 2 . 
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